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ABSTRACT 


1.0. Introduction 


Time-lapse seismic monitoring of production-induced changes in a reservoir and the surrounding rocks over time 
has the basic aim of mapping reservoir compartments and subsurface rock deformation, identifying by-passed oil 
and pore pressure changes, monitoring fluid movement and planning for future production performance (Calvert et 


al., 2018; Herwanger and Koutsabeloulis, 2011; Landrg, 2001; Trani et al., 2011; MacBeth et al., 2018). 


There are the technology, business and geohazard drives that may motivate a research on the effect of hydrocarbon 
production on surface, subsurface structures and faults reactivation using seismic time-lapse. Reservoir compaction 
improves production performance (Bruno, 2002; Setarri, 2002). This is possible because compaction causes 
reduction in porosity and permeability which affect production. Porosity loss affects the computation of oil and gas 
reserves (Pourciau et al., 2005). Using subsidence from time-lapse time shift, subsidence pattern may also be a 
good indication of permeability anisotropy and an indication of bypassed oil (Vedanti et al., 2009). Reservoir 
compaction can lead to reactivation of faults, casing failure and borehole breakout (Morton et al., 2006; Mildren et 


al., 2002). 


Several researchers have made enormous contributions based on time-lapse feasibility studies and reservoir 
monitoring and management within the Niger Delta basin (Aniwetalu et al., 2017; Uko et al., 2018). Moreover, 
literatures in public domain are very scarce on seismic time-lapse, in Nigeria, with the main objective of 
determining effect of hydrocarbon production on subsurface structures and faults reactivation (Igwenagu, et al, 


2021). Ogbamikhumi et al., 2017 evaluated the effect of changes in reservoir fluid saturations on time-lapse seismic 
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amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic 


impeadance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir 
fluid saturation were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where 
porosity is 0.25%. They demonstrated the use of repeat seismic to delineate swept zones and areas hit with water 
override in a producing onshore reservoir. Uko and Otugo, 2016 found based only on orthometric height difference, 
not on reservoir stress changes that the elevation changes in this field are localized only where the measurements 


are located mainly in the river and drenched channels and slopes caused by erosion. 


This study is taken from Fuba Field, Depobelt, Niger Delta, Nigeria. The ultimate deliverable of this study was 
seismic time-lapse effects of hydrocarbon production on the pressure, porosity, and permeability regimes of the 
area. The major components of our study are: (a) Differencing of seismic volume, (b) Well Correlation performed 
in order to determine the continuity of the reservoir sand across the field, (c) Seismic Interpretation which involves 


well-to-seismic ties, fault mapping, horizon mapping, time surface generation and attributes generation. 
“© 2.0. Location and Geology of the Study Area 


The proposed study area Fuba Field is located in the onshore Niger Delta region. Figure 1 shows the map of the 
Niger Delta region showing the study area while Figure 2 shows the base map showing well locations in the study 


area. 
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Figure 1. Map of Niger Delta Showing Figure 2. Base Map Showing Well Locations 
the study Area in the Study Area 


The Niger Delta lies between latitudes 4° N and 6° N and longitudes 3° E and 9° E (Whiteman, 1982). The Delta 
ranks as one of the major oil and gas provinces globally, with an estimated ultimate recovery of 40 billion barrels of 
oil and 40 trillion cubic feet of gas (Adegoke, et al., 2017). The coastal sedimentary basin of Nigeria has been the 
scene of three depositional cycles (Short, and Stauble, 1967). The first began with a marine incursion in the middle 
Cretaceous and was terminated by a mild folding phase in Santonian time. The second included the growth of a 
proto-Niger delta during the Late Cretaceous and ended in a major Paleocene marine transgression. The third cycle, 
from Eocene to Recent, marked the continuous growth of the main Niger delta. A new threefold lithostratigraphic 
subdivision is introduced for the Niger delta subsurface, comprising an upper sandy Benin Formation, an 


intervening unit of alternating sandstone and shale named the Agbada Formation, and a lower shaly Akata 
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Formation. These three units extend across the whole delta and each ranges in age from early Tertiary to Recent. 


They are related to the present outcrops and environments of deposition. A separate member of the Benin 
Formation is recognized in the Port Harcourt area. It is Miocene-Recent in age with a minimum thickness of more 
than 6,000ft (1829m) and made up of continental sands and sandstones (>90%) with few shale intercalations 
(Horsfall. et al., 2017). Subsurface structures are described as resulting from movement under the influence of 
gravity and their distribution is related to growth stages of the delta (Ochoma, et al., 2020). Rollover anticlines in 
front of growth faults form the main objectives of oil exploration, the hydrocarbons being found in sandstone 


reservoirs of the Agbada Formation. 
“= 3.0. Methodology 
3.1. Normalized Root Mean Square (NRMS) Repeatability (RRR) and Differencing of Seismic Volume 


Pre-stack time migrated full-offset 3D and 4D stacks were available. The success of time-lapse reservoir 
monitoring depends on removing the non-repeatable effects such as configurations, seasonal changes, atmospheric 
temperatures, tides, elastic properties of the overburden, compaction, multiples, and rock heterogeneities (Vedanti, 
et al., 2009; Varela, et al., 2006). Obstructions, weather patterns, cost constraints, and maritime traffic can also 
influence the survey orientation. Having considered the above sources of error in repeatability, the acquisition 
system itself, positioning accuracy, receiver sensitivity/calibration, and source calibration must also be looked into. 
All these sources of error were handled through normalized root mean square (NRMS) analysis. The NRMS value 
is simply the RMS amplitude of the difference, normalized by the average of the RMS amplitudes of the base data 
and monitor data (Kragh, and Christie, 2002): 


2 rms (Monitor—baseline) 


NRMS= (1) 
rms (Monitor)+ (Basel ine) 
ye 
ms = ,| = 
af (2) 
The summation is over N number of all samples x;(i = 1, 2, - - -, N) in the time window. 


The NRMS value is a measure of non-repeatability. If NRMS = 0, the data are perfectly repeatable. Typical “good” 
values of NRMS quoted in the literature range from 0.1 to 0.3 [10% to 30% non-repeatability] (Johnston, 2013). 


Monitor 4D seismic volume was subtracted from the Base 3D volume, and the differences was interpreted to 
determine the areas of the field that have been changed during production. Areas of the field where there have been 


changes were analysed and compared to production activity in those areas. 
3.2. Well-to-Seismic Ties 


Well correlation is the first stage of the pre-interpretation process. The process of well correlation involves 
lithologic description, picking top and base of sand-bodies, fluid discrimination and then linking these properties 
from one well to another based on similarity in trends. In between these two lithologies in the subsurface, the 


gamma ray log is often used. Correlation of reservoir sands was achieved using the top and base of reservoir sands 
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picked. The correlation process was possible based on similarity in the behavior of the gamma ray log the Niger 


Delta; the predominant lithologies are sands and shales. In order to discriminate shapes. Also, the thickness of the 
shale bodies overlying and underlying the sand body is considered during Correlation. After defining the 


lithologies, the resistivity log was used for discriminating the type of fluid occurring within the pores in the rocks. 


There are seven basic steps involved in seismic interpretation relevant to this study and they include; Differencing 
of seismic volume, Well-to-seismic ties, Fault Mapping, Horizon mapping, Time surface generation, Depth 
Conversion and attributes generation. Well-to-seismic tie is a process that enables the visualization of well 
information on seismic data. For this process to be achieved, the following are basic requirements; checkshot, sonic 
log, density log and a wavelet. The sonic log, which is the reciprocal of velocity, was calibrated using the checkshot 
data. The calibration process is necessary in order to improve the quality of the sonic log because the sonic log is 
prone to washouts and other wellbore related issues. The results of calibrating the sonic log with the checkshot 


gives a new log called the calibrated sonic log. 


The calibrated sonic log is used along with the density log to generate an acoustic impedance (AJ) log. The acoustic 
impedance log is calculated for each layer of rock. The next step involves generating the reflectivity coefficient 
(RC) log. The RC is calculated and generated using the AI log. The RC log generated is then convolved with a 
wavelet to generate a synthetic seismogram which is comparable with the seismic data. The statistical wavelet 
utilized for convolution is extracted from the seismic data. The synthetic seismogram was generated for every well 
that had checkshot, density and sonic log. The reflections on the synthetic seismogram were matched with the 
reflections on seismic data. The mathematical expressions that governs the entire well-to-Seismic tie workflow are 


presented below; 


Al = pv (3) 
RC = P2V2—Piy4 (4) 
P2V2t+P4V4 
. : P2V2—P1V4 
Synthetic Seismogram = ————— * wavelet 5 
» 2 P2V2t Pi, (6) 


Where p= Density, v =Velocity, AI= Acoustic impedance and RC= Reflection coefficient 


Faults were identified as discontinuities or breaks in the seismic reflections. Faults were mapped on both inline and 
cross-line directions. Horizons are continuous lateral reflection events that are truncated by fault lines. The horizon 
interpretation process was conducted along both inline and crossline direction. At the end of the horizon mapping, 
a seed grid is generated which serves as an input for time surface generation. Time surfaces were generated using 
the seed grids gotten from the horizon mapping process. The third order polynomial velocity model was generated 


and used to depth convert the time surfaces of the reservoirs of interest. 
3.3. Variance (Edge Detection) Method 


The variance attribute is edge imaging and detection techniques. It is used for imaging discontinuity related to 
faulting or stratigraphy in seismic data. Variance attribute is proven to help in imaging of channels, fault zones, 


fractures, unconformities and the major sequence boundaries (Pigott, et al, 2013). In the Petrel software, the 
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variance attribute uses an algorithm that computes the local variance of the seismic data through a multi-trace 


window with user-defined size. The local variance is computed from horizontal sub-slices for each voxel. A 
vertical window was used for smoothing the computed variance and the observed amplitude normalized. The 
variance attribute measures the horizontal continuity of the amplitude that is the amplitude difference of the 


individual traces from their mean value within a gliding CMP window. 


gee y (x, +x,,)° (6) 


Nn f=1 


Where o =standard deviation, o =variance, m =the number of observations, f, = frequency, x, =the variable 


and x,, =mean of x, 


3.4. Dip magnitude 


Dip magnitude and is analogous to strike and dip of sedimentary layers. The dip magnitude is defined as the angle 
between the steepest direction of a plane and a horizontal plane, where values range from 0 to 90. The dip 


magnitude attribute computation in Petrel software makes use of the inbuilt formula: 


tan( 0, ) 


True dip = tan’ (———+ 
tan(f.) 


(7) 


where @, = apparent dip in a direction (y) and #, = dip azimuth relative to a direction (x). 


3.5. Determination of Relative Acoustic Impedance 


The acoustic impedance inversion transforms the seismic data into an acoustic impedance model. The acoustic 
impedance of a media is given as equation (3). To measure acoustic impedance, it is necessary to use seismic 
inversion. It was assumed that the input seismic data has been processed to reduced noise and multiples, and also 


contains zero phase and large bandwidth. The seismic trace represents a band-limited reflective series; 


1 
f@O= 4° (8) 
2 pv 


where f(t) =seismic trace, aOU=the difference in the product of density and velocity. The integration of the 


seismic trace will provide a bandlimited estimate of the natural log of the acoustic impedance. Since the integration 
of band-limited, the impedance will not have absolute magnitude values and consequently is only relative. Relative 
acoustic impedance is an estimated inversion computed by the integration of seismic trace accompany by a high cut 
Butterworth zero-phase filter. It is a simplified inversion and has been generated as an asynchronous attribute in the 
software. It enhances acoustic impedance contrast boundaries. The relative acoustic impedance (RAJ) can be 
computed by integrating the real part of the seismic trace (Taner, 2001). 

t=l 


Inipv)=2f f(Tat ) 


t—0 
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Where f(T) = real part of seismic trace. A Butterworth filter is then applied to remove long-wavelength trends 


that originated from the integration process (Schlumberger, 2007). 


BL( f)=——— 
1 2N 
— 


where BL(f) = band —limited signal in frequency; /,, = frequency cut-off value of 10 Hz, N =filter order of 3. 


(10) 


It is used for delineating sequence boundaries, unconformity surfaces, and discontinuities. The acoustic impedance 


may be related to the formation porosity and the presence of fluid in a hydrocarbon reservoir. 
== 4.0. Results and Discussion 
4.1. Production Data 


The production and reservoir pressure reports are presented in Figures 3 and 4. Production decline of 


1192.21bbl/day resulted from pressure decline of 95.50bar/year or 0.062psi/year. 
4.2. Normalized Root Mean Square (NRMS) Repeatability (RRR) and 4D Response 


The NRMS of 0.38 has been achieved, in this study, implying very good repeatability when considering the quality 
of data and acquisition difference. The seismic time-lapse difference between the base and monitor surveys was 
successfully extracted. The fact of the difference from the monitor implies the existence of production induced 


effect and acquisition, environmental and processing noises, hence the 4D or time-lapse response signal (Figure 5). 
4.3. Reservoir Identification, Correlation and Well-to-Seismic Ties 


The results for lithology and reservoir identification are presented in Figure 6. A total of four sand bodies (L, M, N 
and O) were identified and correlated across all seven wells in the field. Three reservoir sands were selected for the 
purpose of this study (Reservoirs M, N and O). The resistivity logs which reveals the presence of hydrocarbons 
were used to identify the hydrocarbon bearing sands. On Figure 6, the sands are coloured yellow while shales are 
grey in colour. The results for well-to-seismic tie conducted on Fuba field using density log, sonic log and 
checkshot of Well-1 is presented in Figure 7. An extended white 2 wavelet was used to give a near perfect match 


between the seismic and synthetic seismogram. 
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Figure 3. Reservoir Production Rate Versus Days of Figure 4. Effect of Reservoir Pressure Decline on 
Oil Production Production 
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Figure 5. Base, Monitor and Difference Seismic for Figure 6. Well Section Showing Reservoir Section 


Inline 8590 Interpreted Identified and Correlated Across Fuba Field 
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Figure 7. Synthetic Seismogram Generation and Figure 8. 4D Effect (Diffference) Betweeen Base and 

Well-to-seismic tie Conducted for Fuba Field using = Monitor Interpreted Faults Displayed on the Variance 
Well-1 Checkshot Time Slice 


Figure 9. Third Order Polynomial Velocity model for 


Figure 10. Reservoir Surface for Depth Surface M 


Converting Reservoir Surfaces from Time to Depth 


ISSN: 2582-3981 [86] OPEN 6} ACCESS 


& 


Ss Irish Interdisciplinary Journal of Science & Research (IIJSR) 
Ii Js RK Volume 8, Issue 1, Pages 80-92, January-March 2024 


Elevation depth (m] 
10800.00 


Seale. 
178125 


Contour ine 
Ey 


Username 
we 


. re gegeeee 
Br ign 
~S | “4 
sae 
ear ange 
+t » 


Pecans § 
“er e) 


Figure 13. Dip Magnitude Inline 8515 for Base, Figure 14. Variance Edge Inline 8515 for Base, 


Monitor and Difference Volume Monitor and Difference Volume 


Figure 15. Relative Acoustic Impedance Inline 8515 for Base, Monitor 


and Difference Volume 
4.4, Fault Interpretation and Structural Patterns 


The results for the interpreted faults in Fuba field are presented in Figures 5 and 8. Figure 5 shows both synthetic 


and antithetic faults interpreted along seismic inlines for the base and monitor seismic. Faults are more visible 
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along the inline direction because this direction reveals the true dip position of geologic structures. Figure 5a shows 


the base interpreted faults, Figure 5b shows the monitor interpreted faults while Figure 5c shows the difference 
volume with the areas in red circles which show the 4D response. The variance time slice was used to validate the 
interpreted faults as seen on Figure 8. All interpreted faults are normal synthetic and antithetic faults. A total of 
thirty-six faults were interpreted across the entire base seismic data while a total of forty faults were interpreted 
across the entire monitor seismic data. Of the thirty-six interpreted faults, only F4 (synthetic fault) and F13 
(antithetic fault) faults are regional, running from the top to bottom across the field for the base while only F1 
(synthetic fault) and F4 (antithetic fault) faults are regional, running from the top to bottom across the field for the 
monitor seismic. Hence, these faults play significant roles in trap formation at the upper, middle and lower sections 
of the field. There is depletion of hydrostatic pressure due to production. So, areas that were cracks or fissures on 
the base seismic became relaxed and became well-defined or pronounced faults on the monitor seismic. In Fuba 
Field, the fault lengths, dips and orientations are unequal in both base and monitor seismic sections suggestive of 


fault reactivation from hydrocarbon production. 


The results for the interpreted seismic horizons (Horizon M, Horizon N and Horizon O) are also presented in Figure 
5. On these horizons, the fault polygons were generated and eliminated. The horizons were used as inputs for the 
generation of reservoir time surfaces. The reservoir time surfaces (M, N and O reservoirs) reveal that the reservoir 
structure is a collapsed crest, bounded by two regional faults (F1 and F4). Reservoir M, N and O time surfaces are 
truncated by two bounding faults and three minor inter-reservoir faults. The fault supported collapsed crestal 
structures are the closures identified on reservoir M, N and O reservoir surfaces respectively. The similarity in 
structure on reservoir M, N and O reveals that the field is structurally controlled by faults. The third order 
polynomial velocity function used for depth converting the studied surfaces is presented in Figure 9. The depth 
converted reservoir M, N and O surfaces are presented in Figures 10, 11 and 12 for the third order polynomial 
velocity function. The depth structure maps reveal that the reservoirs are anticlinal and fault supported. Reservoir 
M is found at a shallower depth from 10937 to 10997 ft, reservoir N is found at a depth ranging from 11213 to 
11241 ft while reservoir O is found at a deeper depth ranging from 11681 to 11871 ft respectively. 


4.5. Seismic Attributes 


Figure 13 shows the dip magnitude of the faults for the base, monitor and difference volume. The dip magnitude 
values range from 0 to 90 degrees. On Figure 13, the green colors represent areas of greater dip, while the red 
colours represent areas of shallower dip. The dip magnitude of the faults are not equal as indicated by the difference 


volume (Figure 13c). This indicates that there is fault re-activation in Fuba field due to production. 


Figure 14 shows the computed variance attribute of the seismic section. The range of the variance values is from 0.0 
to 1.0. Discontinuities are represented by values of variance equal to 1 while the value of 0 represents a continuous 
seismic event. Red to yellow colourations denote the high values. There are more discontinuities in Figure 14c 
which implies that there are more cracks in the field of study due to production. On the variance map, the areas 
dotted with blue, green, orange and pink colored lines signify values that correspond to the location of the 


discontinuity. The discontinuities may be interpreted as faults and boundaries as shown by the lines drawn on both 
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attribute maps (Law and Chung, 2006). The variance edge enhanced the faults or sedimentological bodies or faults 


within the seismic data volume. Futhermore, high reflectivity sediments and the possibility that a potential 
hydrocarbon trap exists in the area are indicated by several bright spots delineated (in black circles and black ovals). 
Channels, fault zones, fractures, unconformities and the major sequence boundaries can be imaged using the 
variance attribute since it is edge imaging and detection technique. Due to similar seismic traces, there are areas 
with low variance values. Lineaments/discontinuities are indicated by areas with red patches while the structural 


framework of the field is indicated by grey areas. 


The relative acoustic impedance generated in the study area is shown in Figure 15. Based on the map, the highest 
relative acoustic impedance is indicated by the yellow and red colours (in black circles) while the blue colour 
indicates the lowest relative acoustic impedance. Apparent acoustic impedance or physical property contrasts is 
represented by the relative acoustic impedance attribute. It is usually used for sand/shale discrimination, thickness 
variation and indicating sequences boundaries associated with high contrasts in acoustic impedance values, 
unconformity surfaces, discontinuities, porosity and hydrocarbon-bearing sand units (Oyeyemi and Aizebeokhai, 
2015). Shalier facies correspond to high relative acoustic impedance values and may be interpreted as sequences 
boundaries while sand intervals correspond to low values (Alabi and Enikanselu, 2019). Figure 15c shows the 
lowest relative acoustic impedance which indicates a reduction in porosity and permeability due to hydrocarbon 


production. 
“= 5.0. Conclusion 


Production decline of 1192.21bbl/day resulted from pressure decline of 95.50bar/year or 0.062psi/year. Most of the 
data points were not repeatable as evidenced in the computed Normalized root mean square (NRMS) of 0.38 
meaning that on 62.0% of Base and Monitor data points were coincident. A total of four sand bodies (L, M, N and 
O) were identified and correlated across all seven wells in the field. Three horizons (M, N and O) were selected for 
the study. Reservoir M is found at a shallower depth from 10937 to 10997 ft, reservoir N is found at a depth ranging 
from 11213 to 11241 ft while reservoir O is found at a deeper depth ranging from 11681 to 11871 ft respectively. 
Structural interpretation of seismic data revealed that the field is highly faulted with synthetic and antithetic faults 
which are in line with faults trends identified in the Niger Delta. Fault and horizon interpretation revealed that 
closures found are collapsed crestal structures bounded by two major faults. The depth structure maps reveal 
anticlinal faults. The synthetic and antithetic faults act as good traps for the hydrocarbon accumulation in the study 
area. The variance edge enhanced the faults or sedimentological bodies within the seismic data volume. There are 
more discontinuities in the difference volume variance edge which implies that there are more cracks in the field of 
study due to production. The lengths, dips and orientations of the faults and horizons, in the base and monitor 
stacks, are not equal indicative of faults reactivation that could have resulted from hydrocarbon production. The 
difference volume relative acoustic impedance shows the lowest relative acoustic impedance which indicates a 
reduction in porosity and permeability due to hydrocarbon production. In reservoirs, hydrocarbons were 
encountered by all seven wells drilled in the field. Deterministic hydrocarbon volume estimation of the area of 


study indicates large quantities of hydrocarbon at the intervals where these structures were mapped. Within the 
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limits of the available data, it is recommended that further studies should include integration of stratigraphic data of 


all the wells. This will provide more reliable data for interpretation of the depositional environments. Also, further 
studies are required to confirm if these faults are sealing or leaking in which case they serve as conduits for 


hydrocarbon migrations rather than lateral barriers to hydrocarbon escape. 


ISSN: 2582-3981 OPEN @ACCESS 


& 


J S R Trish Interdisciplinary Journal of Science & Research (IIJSR) 
(a oh Volume 8, Issue 1, Pages 80-92, January-March 2024 


SS THROUGH RESEARCH 


ISSN: 2582-3981 OPEN @ACCESS 


& 


J Ss R Trish Interdisciplinary Journal of Science & Research (IIJSR) 
_— a a nk Volume 8, Issue 1, Pages 80-92, January-March 2024 


ISSN: 2582-3981 [92] OPEN@ACCESS 


